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Abstract 

In this study, an algorithm for computing the inverse of periodic k 
banded matrices , which are needed for solving the differential 
equations by using the finite differences, the solution of partial differential 
equations and the solution of boundary value problems is obtained and the 
inverses of periodic anti k banded matrices are computed. In addition, the 
determinant of these type of matrices and the solution of linear systems 
having these coefficient matrices are investigated. When obtaining this 
algorithm, the LU factorization is used. The algorithm is implementable 
to the CAS (Computer Algebra Systems) such as Maple and Mathemat- 
ica. 
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1 Introduction 

The nxn periodic k banded and anti k banded matrices take the following form 
respectively: 
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The periodic k banded matrices are needed in science and engineering 
applications for example solving differential equations by using the finite 
differences, the solution of partial differential equations and the solution of 
boundary value problems. 

In [T], the authotrs obtained an algorithm to find the inverse of the 
periodic tridiagonal matrix by using Doolittle LU factorization and the 
inverse of periodic anti-tridiagonal matrix is obtained when the inverse of 
periodic tridiagonal matrix exists. A new algorithm is obtained for the inverse 
of periodic pentadiagonal and anti-pentadiagonal matrix in [2]. This paper 
is an expansion of [1,. In ^ an algorithm for solving linear systems having 
periodic pentadiagonal coefficient matrices is obtained. It is presented that a 
new computational algorithm to evaluate the determinant of the tridiagonal 
matrix with its cost In [4 . In [5] Hadj and Elouafi obtained a fast numerical 
algorithm for the inverse of a tridiagonal and pentadiagonal matrix. 

In this work we obtain an algorithm to find the inverse of periodic k banded 
matrix when its inversion exists. When the algorithm is obtained the Doolittle 
LU factorization is used. After finding the inverse of periodic k banded matrix, 
the periodic anti k banded matrix is inverted by using the inversion of periodic k 
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banded matrix. Also, an algorithm is studied to solve the linear systems having 
these coefficient matrices. 



2 Main Result 

In this section, the LU factorization of the matrix M is computed firstly where 
L and U are lower and upper triangular matrices, respectively. It is as in the 
following: 
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for r = i — 1, i — 2, . . . i — and 
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for i,i < then lij = and Uij = 0. 
From here, it is clear that 

n 

det{M) =]Juii. 

i=l 

The inverse matrix is computed as in the following, if the matrix M is 
nonsingular: 

Let CrS be the rth column of for r = 1,2, . . . ,n then 

M ^ — {Sij)i<ij<:n = (Cl, C27 • • • 7 C'ri • • • I Gn).Cr = ('S'l,r5 '5'2,r) • • • ) "^nir) 

and it can be written as in the following: 

Cr = (Ci , C2, . . . , Cr, . . . , Cn)Er 

where is the Kronecker symbol {Ej. = ((5ir,^2r7 • • • 7 ^nr.)'^ , r ~ 1,2, . . . , n.). 

Now, the algortihm for the inverse of the periodic k banded matrix can be 
developed. By using the LU factorization, the entries of the last columns 
of are computed as follows: 
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,^ = r-l,r-2,...,2,l 



where ti^s are the entries of the inverse of L and for r = l,2,...,n — 1 they are 
computed with the following reccurence relation: 
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Up to now, the entries of the last -^^^ columns are obtained. Entries of the 
2 



remaining n — columns are computed by using the following equation 
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where if j = 1,2, — then a-^^k^ / 0. Here, i > n, j > n then 

= 0, if i > n then Ci — 0. 
Algorithm 1: 

INPUT: n is the order of the periodic k banded matrix, k is the bandwidth 

of the matrix. 

OUTPUT: The inverse matrix Af-^ = (S'jjOk, j<„ • 
Stepl: For i = 1, 2, . . . , — ^ii, if = 0, then i^h^ ~ ^■ 

Step2: For i = ... ,n, if a- ^^k^ = 0, then a- A. 

Step3: By using (O-©, compute the elements of the matrices L and U. 

For i = 1,2, ... ,n, a Uii ^ 0, then ua = X. 

Step 4: Compute det(M) = ^ Uii . If M is singular, then the output 

\i=l / A=0 

is " Singular Matrix" . 

Step 5: Compute the elements ti,r by using 

k+l 



Step6: Compute the elements of the last columns by using ([8]). 
Step7: Compute the elements of the remaining n — columns by using 
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StepS: Substitute the actual value of A in all elements of the inverse matrix 

M-\ 

The inverse matrix of the periodic anti k banded matrix N can be obtained 
by using the inverse of periodic k banded matrix M easily. 
Let i? be an n X n matrix as in the following form: 
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It is clear that R is nonsingular and its inversion is itself [2]. 
The following relation is true for (??)and 

N = MR 

Thus, the inverse of © is obtained as in the following: 
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Also, linear systems having periodic k banded coefficient matrix, Mx = y, 
can be solved by using the Dolittle LU factorization of this type of matrix. Here 
the component uu is important. Because, the system has a unique solution if 
Uii 0, i = 1, 2, . . . , n. 

Algorithm 2: 

Stepl: By using the Stepl-Step3 of Algorithml, compute the LU 
factorization of the coefficient matrix M. 



Step2: For i — 1,2, . . . ,n Zi — yi 
StepS: For i — n,n ~ I, ... ,2,1 
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3 Numerical Example 

Example 1 Consider the 6x6 matrices M and N as in the following 
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We apply the Algorithml to the matrix M and we have 

• For i = 1, 2, 3, 4 = (^12, W23, W34, W45) = (1, 2, 3, 1) 

• Fori = 2, 3, 4, 5 = (/21, ^32, ^43, ^54) = (5, 
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Ci = {Sei, S51, 84,1, S31, S21, Sii) = (— 1,0, ■ 
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Example 2 Consider the 10 x 10 matrices M and N as in the following 



M 



N = 



1 


-1 


2 


2 


-1 














1 


2 


-1 


3 


1 


1 


2 














1 


-1 


1 


2 


1 


-2 


-1 











3 


1 


-1 


1 


-3 


1 


1 


-3 








2 


-1 


1 





-3 


2 


1 


-1 


-1 








1 


2 





-1 





-2 


1 





1 








-2 





1 


-1 


1 


-2 


1 


-1 











1 


3 


2 


-1 


1 


2 


1 














-1 





2 


1 


-2 


1 


2 














2 


1 


1 


-1 


2 


2 














2 


1 


1 


-1 


2 














-1 





2 


1 


-2 


1 











1 


3 


2 


-1 


1 


2 


1 








-2 





1 


-1 


1 


-2 


1 


-1 





1 


2 





-1 





-2 


1 





1 


2 


-1 


1 





-3 


2 


1 


-1 


-1 





3 


1 


-1 


1 


-3 


1 


1 


-3 








1 


-1 


1 


2 


1 


-2 


-1 











2 


-1 


3 


1 


1 


2 














1 


-1 


2 


2 


-1 














1 



We apply the Algorithml to the matrix M and we have 
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• Fori = 1,2,3,4,5,6,7,8 

Ui,i+i = (^12, U23, U34, U45, useuer, urs, ugg) = (-1, -1, 0, 6, 3, -7, -||, ^) 

• Fori = 1,2,3,4,5,6,7 
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Example 3 Solve the periodic pentadiagonal system given as follows: 
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We apply the Algorithm2: 

• Fori =1,2, 3, 4 = (ui2, U23, "34, M45) = (2, -1, |, -7) 

• For i = 1, 2, 3 = (wis, "24, W35) = (-1, 1, 2) 
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• For i = 2, 3, 4, 5 h^i-i = {hi, hi, ks, hi) = (2, 5, 3, -2) 

• Fori = 3, 4, 5 Zi,i-2 = («3i, ^42, ^53) = (1, -|, -5) 

• For i = 1,2,3,4,5,6 Ujj = (un, U22, W33, W44, U55, ^ee) = (1, -5, |, -1, -3, 

• For i = 1, 2, 3, 4, 5 Uj,6 = (wi6, "26, W36, W46, wse) = (1, -2, -f , -1, -2) 

• For z = 1, 2, 3, 4, 5 Ze.i = {hi,h2, ha, hi, hb) = (1, |, 7, 5, -f ) 

• det(M) = 14 

• {zi,Z2,Z3, Z4, Z5, zb) = (3, -7, ^, -9, -5, -|) 

• (a;i,a;2,a;3,X4,X5,X6) = (1, 1, 1, 1, 1, 1) 
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